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Convergence of approximate solution 
of delay Volterra integral equations 


M. Zarebnia * and L. Shiri 


Abstract 


In this paper, sinc-collocation method is discussed to solve Volterra func- 
tional integral equations with delay function 0(t). Also the existence and 
uniqueness of numerical solutions for these equations are provided. This 
method improves conventional results and achieves exponential convergence. 
Numerical results are included to confirm the efficiency and accuracy of the 
method. 


Keywords: Volterra functional integral equations; delay function; sinc- 
collocation. 


1 Introduction 


Delay integral equations arise widely in scientific fields such as physics, 
biology, ecology, control theory, etc. Due to the practical application of 
these equations, they must be solved successfully with efficient numerical ap- 
proaches. In recent years, there have been extensive studies in convergence 
properties and stability analyses of these numerical methods, see, for exam- 
ple, [10]. The numerical solutions of integral equations with delays have also 
been discussed by several authors such as Brunner [1], Li and Kuang [5], Linz 
and Wang [6]. 

Sinc methods for approximating the solutions of Volterra integral equa- 
tions have received considerable attention mainly due to their high accuracy. 
These approximations converge rapidly to the exact solutions as the number 
of sinc points increases. Systematic introduction of these methods can be 
found in [9]. In [11] sinc-collocation method is emplyed to solve Volterra 
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integral equations with smooth kernels. The analytical and numerical tech- 

niques used in these works can be extended to delay integral equations. 
The main objective of the current study is to implement the sinc- 

collocation method for Volterra functional integral equation of the form 


y(t) = g(t) + (Vy)(t) + Voy) (4); teT:= [0,7]. (1) 


The Volterra delay integral operators V and Vg (from C(I) + C(Z)) describ- 
ing these equations are defined by 


(y(t) == i Ky (t,)y(s)ds 


and 


6(t) 
(Voy) (t) == | Kolt, s)y(s)ds, 


respectively, and the delay function 6 is subject to the following conditions: 
(D1) 6(0) = 0, and @ is strictly increasing on the interval I; 

(D2) 6(t) < qt, t € I, for some ¢ € (0,1); 

(D3) 6 € C4(I) for some d > 0. 

We will refer to a @ that satisfies (D1) as a vanishing delay function (or, in 
short, a vanishing delay). The linear case, 0(t) = gt =t—(1—q)t =:t—7(t) 
(0 < q < 1) (proportional delay) is also known as the pantograph delay 
function [4]. In this paper we consider vanishing delay but our methods can 
be use with nonvanishing delay too. 

The layout of this paper is as follows. In Section 2, the solvability of 
equation (1) is stated. Section 3 outlines some of the main properties of sinc 
function that is necessary for the formulation of the delay integral equation. 
Sinc-collocation method is considered in Section 4. In section 5, we analyze 
the existence and uniqueness of numerical solutions. In Section 6, the order 
of scheme convergence using the new approach is described. Finally, Section 
7 contains the numerical experiments. 


2 Existence and uniqueness of solutions 


In the present section, we state the solvability of integral equations with 
vanishing delay. The following theorem generalizes Volterras 1897 classical 
result on the existence and uniqueness of solutions for the equation (1) with 
O(t) = qt (0<q< 1). 


Theorem 1. /3] Assume that the given functions g, Ky, and Kg in (1) satisfy 
1) g€ C(I), Ki € C(D), and Kz € C(Do), where 
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D:={(t,s):0<s<t<T}, Do:={(t,s):0<s< 6),teT}; 


2) @(t) is subject to the assumptions (D1)-(D3). 
Then for each g € C(I) there exists a unique function y € C(I) which solves 
the equation (1) on I. 


3 Review of the sinc approximation 


In this section, we will review sinc function properties, sinc quadrature rule, 
and the sinc method. These are discussed thoroughly in [9]. For any h > 0, 
the sinc basis functions are given by 


z—gjh 


S(j, h)(z) = sinc( 


ds j =0, 1, Diss atey 


where 


sin(7z) ‘ 
sinc(z) = { mz L aa ie 
ie OP 


The sinc function form for the interpolating point z, = kh is given by 


sunyeen) ={ 9p o9 


They are based in the infinite strip Dg in the complex plane 

Da = {w=utiv: |r| < dh. 
To construct approximation on the interval [0,7], we consider the conformal 
map 


z& 


o(2) = In 


). 
The map ¢ carries the eye-shaped region 


Zz 
D= : — d}. 
{z€C: |arg( rape } 
The function 
Te” 
Soy d. = 
z=¢ (w)= os 


is an inverse mapping of w = ¢(z). We define the range of ¢~! on the real 
line as 


T= {¢(u) = 6 "(u) € D: -00 <u < oo}. 
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The sinc grid points z, € (0,7) in D will be denoted by x, because they 
are real. For the evenly spaced nodes {kh}? ., on the real line, the image 
which corresponds to these nodes is denoted by 


Tekh 


tp = 6 (kh) = a k= +1, 


Definition 1. Let D be a simply connected domain which satisfies (a,b) C D 
and a and c; be a positive constant. Then £_4(D) denotes the family of all 
functions u € Hol(D) which satisfy 


|u(z)| < e|Q(z)|* (2) 


for all z in D where Q(z) = (z — a)(b— 2). 
The next theorem shows the exponential convergence of the sinc approx- 
imation. 


Theorem 2. Let u € £La(D), let N be a positive integer, and let h be selected 


by the formula h = 
of N, such that 


ae then there exists positive constant cz, independent 


sup |u(t) ee (d(t))| < coVNe~V™ON 


tE (a,b) 
The error analysis of the sinc indefinite integration has been given in [7]. 


Theorem 3. Let uQ € La(D) ford withO<d<7m. Leth = \/24. Then 


there exists a constant Co, which is independent of N, such that 


sup fu nde YS oe, h)(dsw(t))] < ese" VN), (8) 


t€(a,b) 


where 


*—J sin(rt) 


I, h)(x) = 5 rf a dt. 


4 Sinc-collocation method 


In this section, we apply sinc-collocation method to solve equation (1) which 
we state again for the convenience of the reader: 


6(t) 
@+f Ky(t, s)y(s)ds + Ko(t, s)y(s)ds 
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if t = 0 we have y(0) = g(0). For ease of calculation, we employ the trans- 
formation 


T-t 
t) = y(t) - ——- 
u(t) = y(t) — 909), 
in this case u(0) = 0. Then the above problem becomes 
6(t) 
@+f Ky(t, s)u(s)ds + Ko(t, s)u(s)ds (4) 


where 


F(t) = g(t) 


1 1 i 9(t) 
ae t)g(0) 4 790) i, Ki(t,s)(T — s)ds+ : Ko(t, s)(T — s)ds >. 


Now, let u(x) be the exact solution of (4) that is approximated by the 
following expansion 


N 
oh h)oo(t) + u(tw+1) w(t), (5) 


we choose w(t) so that above formula interpolate function u at the points t,, 
so 


N 
1 
w(t)h= alt D7 4 SG,4)(@@) 
j=—N 
where the points t; are defined by 
of Tj ain a i, 


By replacing approximate solution (5) in t = t, in the equation (4), it follows 
that 


N 
S- ujS(j,h)(P(te)) + ungiw(te) 
=—N 
N tg ‘ 

= [PO Kults.s)8.1N(6(9))ds Fuses [Kult s)o(s)as 
N 
ye 


O(tk) (tx) 
uf Ko(tz, s)S(, hy(o(s))ds-+ unr f Ko(tr, s)w(s)ds 


+ flts). (6) 
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We are interested in approximating the integral in above equation by the 
quadrature formula presented in (3). Then by using Theorem 3, we obtain 


FC GEN os. 2 
Bey 7 Molt) k=-N,...,N 41, 


From definition of t; we can write 


J(j,h) (Kh), 


J Kilts, s)SU.Woo(s)ds ~ 
0 


(i AY(6lte)) = { 
The analogue of above equation we have 
(tr) 
J Kaltss)S(G.h)oo(s)ds = 
0 


in which ¢; := $(0(t,)), in next section these formula will be discussed. 
Finally, let 


tk 
Aik — | Ki (tr, s)w(s)ds, 
0 


(te) 
Koz =} Ko(t,, 8)w(s)ds, 
0 
by = Kip + Kor, k=—-N,...,N+1. (7) 


By using relation (3), we can approximate b; in the following form 


th 
Kir = | Ki (tg, s)w(s)ds 
0 


1 tk N 
ag [ Kilt, 8) s— 2, 5G. MOB) ds 
1 tr N th 
= : shag D8 i Ki(t, 8) S(j,h)(4(s))ds 
1 tk N 1 
Be [ eta Gade te ia eat) gay Oe) (8) 


Thus equation (6) is written as 


>? ae) yA testy) IG, h)(Olte)) + Bolte, ty) IG, h) (br) } Uy 


eos = (9) 


Convergence of approximate solution of delay Volterra integral equations 45 


This linear system of equations is equivalent to (4). By solving this system, 
the unknown coefficients u; are determined. We rewirte the linear system 
(7) in matrix form 


[Z-AJUy =F (10) 
where 


An = {41 (te, tj) J, h)(O(te)) + Kaltes ty) IG, bh) (Ge) } 


k=—-N,...N+1, j=—N,...N, 


Ga i) 


Ax,w4i = dr, k=-—N,...N +1, 
Vitae: P= VOR SI 


5 Existence and uniqueness of the sinc-collocation 
solution 


In this section, we study the existence and uniqueness of the solution to (8). 


Lemma 1. For x ER, the function J(j,h)(x) is bounded by 
|G, h)(@)| < 1A. 


Theorem 4. Assume that K,, Ko and f in the Volterra integral equation 
(4) are continuous on their respective domains D, Dg and I. Then there 
exists anh > 0 so that for any h € (0,h) the linear algebraic system (8) has 
a unique solution Un. 


Proof. We know that 


Ang = {41 (te, t5) J, h)(O(te)) + Kaltes ty) IG, h) (Ge) } 


WN Jo Nae. 2g Ns SNe 
Az,w4i = dp, k=—N,...N +1, 


ve i) 


so we can write 


l|Alloo = ,_ max, fh 2 ag 5 | 1s) Ol te)) Ka (tasty) + Kaltes ty) IG A)(Ge)| + IPal}- 


Using Lemma 1, we have 
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oo S Ki (tr, Ka(t, by| }- 
Allo < om ny aa) | 1 (te, tj) + Ka(te, tj)| + bel } 
Theorem 3 and continuty Ky and K2 and equations (7) and (8) give 


bx < ce ¥ TaN | 


Therefore 
|| Aloo < 11h » Fe) 5 |Kultest t;) + Ka(th,t;)| + ce V7, 


Thus the elements of the matrice A are all bounded. The Neumann Lemma 
then shows that the inverse of the matrix J — A exists whenever ||.Al].. < 1. 
This clearly holds whenever h is sufficiently small. In other words, there is 
an h > 0 so that for any h < h matrix A has a uniformly bounded inverse. 
The assertion of Theorem 4 now follows. 


6 Convergence analysis 


The convergence of the sinc-collocation method which is introduced in the 
previous sections is discussed in the present section. It is assumed that wu is 
the exact solution of Eq. (4) and U/y is an approximation of the sinc method. 
Firstly, we state the following lemma which is used subsequently. 


Lemma 2. ( /8/) Leth > 0. Then it holds that 


N 
sup D> 1SG,h)()| 
wR 5" 


In the following theorem, we will find an upper bound for the error. 


a 


(3+InN). 


Theorem 5. Let Un(x) is the approximate solution of integral equation (4). 
Then there exists a constant cs independent of N such that 


sup |u(x) —Uy(x)| < csVN In Ne~V7™2%, (11) 
xe(0,T) 


Proof. For collocation error e := u—Uy 
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N 


sup |je(t) — S- S(j, h)(O(t))e; + en+1 w(t) 


te(0,T) j=—N 


N 


< sup fe(t)— S> S(j,h)(H(E))ej) + lewaal sup |w(t) 
tE(0,T) ae, te(0,T) 


<VNeV™®N + lev sal= sup |t— > t/S(4,h)(6(8)) 


LT te(0,T) j=—-N 


< co V NeW Vtden + lenuilne ve —vadaN 
< ev Newv 7988 


1. ON 
so we can write 
e(t)= S~ S(j,h)(H(t))ej + ensiw(t) + eV/Ne~V™ON 
j=—N 
For t = t; it satisfies the error equation 


e(t) = (Ve)(t) + (Voe)(t). 


The contribution of V in the above error equation is described by 


Z | eee 


tj 
=) Ky (t;,8 > dX S(j, h)(O(s))ex + en+1w(s )ev evra} 
0 


N 
K t yt —vV17da —V TAM 
=h ss eG h)(b(tk) ex =r Ayjen41 + C2e€ or + tjev Ne eae 
k=—N 


N 
Kilt; A A or 
=h 2s i oe */ I(5,h) (P(te) Jew + Kigensa + do VNeW~V™ON, 
also for delay operator Vg we have 
A(t;) 
oe)(t;) = fo Kaley, s)elsas 
0 


N 
= (t; te —Vrda 
=h d 2( buy : j,h)(bx)ex + Kojeny1 +’ VNe Vide’ 
k= 
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Thus, the representation of e; has the form 


inh Y oi 


+ bjeny1 +ev Ne’ maa 


{Ki (tj, th) J(9,h)(O(te)) + Ka(t;, te) IG, h) (x) } ex 


we may write the collocation equation as 
[Z — Ale = eV Ne V7", 


Here, Z and I denotes the identify matrix and constant vector 1, repectively. 
It thus follows from Theorem 4 that uniform bound exists for (Z—.A)~?, so 


llel| < eV NeW ¥ 740% (12) 


Hence, by Lemma 2 and (12) we can obtain the upper bound (11). 


7 Illustrative examples 


In this section, the theoretical results of the previous sections are used for two 
numerical examples. The numerical experiments are implemented in Matlab. 


Example 1. The pantograph Volterra integral equation 


t 
K(t, s)y(s)ds 
a(t) 


with 0(t) = qt, k(t,s) =ts, and g(t) = (-t? +t + 1)e! + t(qt — 1)e", has the 
exact solution y(t) = e’. The results are shown in Table 1. 


Table 1: Values of ||E||oo for Example 1 


N\q 0.01 0.09 0.1 0.5 0.99 
10 3.7266 x 10-© | 4.5684 x 10-® | 4.1900 x 10-° | 8.8432 x 10-> | 1.7511 x 10-° 
30 5.5006 x 107! | 3.6148 x 10-° | 5.2450 x 10-° | 4.0946 x 10-7 | 7.5530 x 1078 
50 8.6331 x 10718 | 1.0118 x 107-1° | 1.3788 x 10-18 | 1.0852 x 10-8 | 1.1792 x 10-9 
70 7.9936 x 10-' | 5.4152 x 107! | 7.3683 x 107!? | 5.8153 x 107! | 1.0688 x 10719 
90 3.2196 x 10-1 | 4.3609 x 10-18 | 5.9374 x 10-18 | 4.6795 x 10-1! | 8.6006 x 10-12 


Example 2. Consider 
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Figure 1: The errors for g = 0.5 in Example 1 


u(t) = g(t) + | K(t,s)y(s)ds 


100 


with k(t,s) = s—t. We choose g(t) so that its exact solution is y(t) = t—t?. 


Table 2 shows the numerical results. 


Table 2: Values of ||E||oo for Example 2 
N\r 0.01 0.09 0.1 0.5 0.99 
10 1.1098 x 10-© | 4.6373 x 10-5 | 5.3418 x 10-® | 5.2686 x 10-© | 5.5796 x 10-% 
30 4.0592 x 10-1 | 8.1311 x 10-!° | 8.1496 x 10-1° | 7.9947 x 10-19 | 8.1273 x 10710 
50 8.2652 x 10-15 | 1.3631 x 10- | 1.4027 x 10-1? | 1.3289 x 10-! | 1.3898 x 10-12 
70 4.8580 x 10-1 | 7.1871 x 10- | 7.0991 x 10-1 | 6.8972 x 10-! | 7.0499 x 10-15 
90 1.4498 x 10716 | 2.2706 x 10716 | 1.9428 x 10716 | 1.6653 x 10716 | 3.3306 x 10716 


-16 i 


20 


Figure 2: The errors for r = 0.5 
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8 Conclusion 


Several methods has been presented for the numerical solution of equation 
(1) in the special cases for example 0(t) = qt [2]. We propose a numerical 
algorithm in order to solve the delay integral equation (Eq. (1)) where 6 is 
general function. Our method has been shown theoretically and numerically 
to be extremely accurate and achieve exponential convergence with respect 
to N. 
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